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The current-voltage characteristics of the a-Hemolysin protein pore during the passage of single- 
stranded DNA under varying ionic strength, C, are studied experimentally. We observe strong 
blockage of the current, weak super-linear growth of the current as a function of voltage, and a 
minimum of the current as a function of C . These observations are interpreted as the result of the 
ion electrostatic self-energy barrier originating from the large difference in the dielectric constants of 
water and the lipid bilayer. The dependence of DNA capture rate on C also agrees with our model. 
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The voltage-driven translocations of polynucleotides 
through nanoscale pores has been recently studied in 
vitro at the single molecule level 1-, 13] . The dynamics 
of biopolymer translocation through nanopores is central 
to many biological processes such as RNA export and 
phage infection, and is the underlying principal behind a 
number of new methods for nucleic acids analysis. Single 
molecule DNA and RNA translocation experiments, as 
well as theoretical models and simulations have provided 
some critical information on the dynamics of biopoly- 
mer transport and its dependence on physical parame- 
ters such as the electric field intensity, polymer length, 
temperature, and chemical characteristics, such as its se- 
quence |llQ,y|. However, much less is known about the 
nature of the ion current during the translocation of the 
biopolymer through the nanopore. 

We use a single a-Hcmolysin (a-HL) pore embedded in 
an insulating phospholipid membrane. An ion current of 
/ ~ 80 pA is reduced down to Ib ^ 7 pA upon the elec- 
trophoretic threading of single stranded DNA molecule 
into the pore (Fig. 1), with the ratio Ib/I = 0.09 j^. 
Under these conditions the translocation time of a single 
base in the polynucleotide is ~5 /is ^ . To establish a cur- 
rent of ~ 7 pA approximately 220 ions must flow through 
the pore in the opposite direction during the passage of 
a nucleotide. As a first approximation we can therefore 
assume that the DNA is nearly static compared with the 
fast moving ions. 

Why is Ib /I so small? The first possible explanation is 
that the passage of the ssDNA strongly reduces the cross 
section of the channel available for ionic movement. In 
a bulk solution ssDNA has the tendency of base stack- 
ing leading to formation of a helix with a diameter of 
^ 1 nm Q , which is considerably smaller than the av- 
eraged channel diameter 1.7 nm |y]). Therefore, such 
a dramatic current blockage is unlikely to be explained 
solely by the reduced cross section. 

In this paper we study how the blocked current Ib 
depends on ionic strength C and voltage V. We also 
find that the capture rate of DNA in the pore strongly 
depends on C. We explain these observations as a conse- 
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FIG. 1: I-V measurements (a) for an "open" pore, and 
(b) when single stranded DNA is electrophoretically threaded 
through the pore. DNA abruptly blocks the ion current from 
80 pA to ~ 7 pA (Ib), as shown in the inset (C = 1 M KCl, 
V=120mV, r = 8°C). 



quence of an electrostatic self-energy barrier Q related to 
the huge difference between dielectric constant ki = 80 
of water and that of lipids and DNA K2 ~ 2. This differ- 
ence results in the confinement of the electric field lines 
of an ion during the ssDNA passage through the channel, 
leads to a large self-energy of the ion, and electrostatically 
amplifies the effect of the channel narrowing. 

An electrostatic self-energy barrier was predicted by 
Parsegian J3, but its effects have not been directly ob- 
served in the conduction of biological channels. Appar- 
ently, evolution has used several compensating mech- 
anisms to facilitate ionic transport (screening by salt, 
"doping" channels walls by fixed charges and coating of 
walls by hanging dipoles) Q, H, llfl IHI' On the other 
hand, electrostatic barrier was observed in conduction 
of synthetic ion channels [l3|, which do not experience 
evolutionary pressure. Similarly, a-HL is not known to 
function as a DNA transporter in vivo. Thus, one might 
expect to observe effects of the self-energy barrier. 

Current- voltage (/ — V) measurements of a unitary a- 
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FIG. 2; Open pore current characteristics at T = 8°C. The 
main figure displays the dependence of the current on the bulk 
KCl concentration measured at 120mV, ISOmV' and 180mV 
(circles, squares and triangles respectively). Lines are linear 
regression fits. The inset is the measured I-V at 1 M KCl. 



HL were performed using dynamic voltage control [Ij, 
[iJI- Fig- 2 displays the current for the open pore as a 
function of C measured at 120 mV, 150 mV and 180 mV. 
In all cases the current follows a linear dependence on C 
(thin lines), in the measured range (0.25 M - 2M). The 
inset displays the I-V characteristic curve measured at 
C= IM KCl. 

The blocked ion current, Ib (Fig- 3), displays several 
marked differences as compared to the open pore current. 
In contrast to the linear dependence of the open pore 
current on C, the blocked current is not monotonic. It 
grows with C roughly linearly for C > 1 M, while at 
C < 1 M it changes only weakly and goes through shallow 
minimum near C — 0.5 M. Furthermore, the I — V curve 
of Ib measured at IM KCl (inset of Fig. 3) is more non- 
ohmic than that of the open pore current displayed in 
the inset of Fig. 2. 

The dependence of the DNA capture rate on the salt 
concentration is displayed in Fig. 4. Previous studies 
showed that the capture rate scales linearly with the bulk 
DNA concentration under similar conditions Ilq|. The 
capture rate dependence on C is highly non-linear: below 
0.5 M KCl the capture rate sharply decreases, and for 
C > 0.5 M it levels off. This trend is consistent for the 
three different voltages that we measured. Measurements 
below 0.25 M were impractical due to the long delay time 
between events. 

We use a crude model to qualitatively interpret the 
experimental results. Our model is lacking sophistica- 
tion and quantitative approach of modern theories of 
ion channels jsl lHI, but emphasizes two-dimensional 
specifics (see below) of DNA translocation physics, which 
has not been discussed in the literature. 

Let us treat ssDNA and the channel internal wall as 
coaxial cylinders (Fig.[SJ) with radius r = 0.5nm and a~ 



< 
a. 




FIG. 3: Main figure: The dependence of the blocked ion 
current, Ib on bulk ion concentration measured at 120mV , 
150ml/ and 180mV (circles, squares and triangles respec- 
tively). Ib values were determined from the peak of the 
distributions of > 1000 DNA translocation events for each 
C and V. Solid lines are guides to eyes. The inset displays 
the measured I-V curve for the blocked pore, performed using 
dynamic voltage control |13. ll4|l . 
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FIG. 4: The dependence of DNA capture rate (per mole) on 
KCl concentration, C. Circles, squares and triangles corre- 
spond to V = 120mV , 150»nV and 1 80m, V, respectively 



0.85 nm, respectively. Salt ions are located in the water- 
filled space between them, with thickness d = 0.35 nm. 
The length of the channel is L = 5 nm. Each charge on 
a ssDNA phosphate in the channel is typically neutral- 
ized by a K+ ion, forming a pair. In the confined space 
of the ssDNA-occupied pore, K+ ion is strongly bound 
to a phosphate moiety on the ssDNA backbone. The 
strong attraction arises from confinement of an electrical 
field between two cylinders. Within this space opposite 
charges interact with strong two-dimensional logarithmic 
potential (2e^/Ki) \n{p/d) (sec Fig. El, where p>d is the 
two-dimensional distance between them. The fact that 
dc^lB/'^: where Ib = e^ /nikBT = 0.7 nm is the Bjerrum 
length and T is the temperature of the experiment, means 
that T ~ Tc/2, where Tc = e^/{kBKid) is the Kosterhz- 



Thouless transition temperature Il6l|. Thus, the thermal 
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FIG. 5: The side view of the membrane and the channel with 
DNA inserted. A phosphate and a K^ ion bound to it are 
shown. 



motion can not break a pair |17J . 

Considering charge transport we may first ignore neu- 
tral pairs. An electric current is produced by an extra ion 
crossing the channel in the confined space. Such ion has 
a higher self-energy than ions in the bulk solution, or in 
other words it goes through an electrostatic self-energy 
barrier U{x). This blocks the ion current. To estimate 
U{x) we write U{x) = e(j){x)/2, where <j){x) is the electro- 
static potential created by a charge e located at x. Our 
numerical calculation in the limit of infinite ratio ki/k2, 
when all electric lines stay in the channel and at a/d ^ 1 
can be well fitted by the following expression: 




FIG. 6: An unfolded view of the water-filled space containing 
an extra K^ ion. Dashed lines represent the electric field 
fines of the charge. At p < a this electric field spreads in all 
directions and becomes uniform far from the charge. 
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The origin of the two terms in Eq. |^ is illustrated in 
Fig. El for X — 0. At d < p < a the electric field of 
the central charge gradually spreads over all azimuthal 
angles in the the whole water-filled space decaying as 
E — e/{Kipd). This leads to (t){x) = {2e/Kid)\n{a/d) and 
the barrier term U2{x). The U2{x) is essentially inde- 
pendent of X at the distances larger than a from the ends 
and but vanishes at the channel ends, where most of elec- 
tric lines are attracted to the bulk solution (this decay 
is not reflected by Eq. (Q)). On the other hand, f/i(0) is 
created by the one-dimensional uniform electric field at 
distances p>a. For |a::| >0 the electric field at the closer 



end is stronger than that of the other end, therefore Ui{x) 
decreases parabolically with \x\ and vanishes at the chan- 
nel ends [Il|. For L = 5nm, a = 0.85nm and (i = 0.35nm 
Eq. (HJ gives C/(0) ~ A.6kBT, or Ui{0) = 2.9 keT and 
1/2(0) = 1.7 ksT. The total barrier U{x) of an extra K+ 
or Cl^ ion is shown on Fig. [3l by the upper curves. 
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FIG. 7: Energy band diagram for K^ ions (solid lines) and 
Cl~ ions (dashed lines). The lower band represents the energy 
of the cations bound to DNA phosphates. The empty upper 
bands show the self-energy of the extra salt cation (solid line) 
and anions (dashed line) entering the channel, a) In the ab- 
sence of the Donnan layers (C > Cd) b) with Donnan layers 
of the width D creating potential —Ud. Vacant phosphates 
are shown by empty circles. The chemical potential p of K^ 
ions in the system is shown by the thin dotted line. 

Recall that there are K+ ions bound to ssDNA phos- 
phates in the channel. Each of them can be removed to 
the bulk creating a vacancy. The energy penalty for this 
process is the same as the penalty for placing an extra 
ion in the same place. Thus, energies of bound K"*" ions 
are —Ui{x) — U2{x) and can be shown by the lower solid 
curve of Fig. [7^. Vacancies have to overcome the barrier 
C/(0) to cross the channel. 

The large self-energy of extra charges deep in the chan- 
nel results in accurate neutralization of DNA by salt 
cations. Such nearly perfect neutralization was observed 
in molecular dynamics modelling jla] of the channel. 

At small salt concentration C < Cd , some cations close 
to the channel ends can escape to the bulk because the 
cations enjoy larger entropy in the solution. As a result 
there are negative charges in the layer of width D at each 
end, and positive screening charge in the adjacent layers 
of the bulk solution. These double layers of the width D 
(see Fig.[7l3) produce the Donnan potential —Ud in the 



channel, and prevent remaining cations from leaving the 
channel. The Donnan potential moves down energies of 
both cation bands, while bending these bands up in the 
very ends (Fig.[7|D). On the other hand, the energy band 
of an anions is moved up. This leads to the exclusion of 
anions from the channel noticed in Rcf. 18] . 

These ideas can be used to interpret the Ib{C) curves 
in Fig. 131 At large enough salt concentration of the bulk 
solution C > Cd, where Cd '--^ 1 M, ion current through 
the channel should be due to extra salt cations and an- 
ions. Therefore the blocked current is proportional to salt 
concentration Ib ocC. On the other hand, at C <Cd the 
Donnan potential repels anions and the charge transport 
is due to cations only. In the first approximation the cur- 
rent Ib is independent of C, because Donnan potential 
Ud grows with decreasing C, resulting in a smaller bar- 
rier t/(0) and for cations to compensate for the decreasing 
C. In the second approximation, at C<C'd the Donnan 
potential Ud changes slightly slower than the chemical 
potential fi = kBT\n{C/CD)- As a result when C de- 
creases the barrier for vacancies decreases slightly and the 
barrier for extra cations slightly increases. This leads to 
vacancy dominated transport [19j and explains the weak 
increase of Ib {C) as C decreases shown in Fig. O In the 
limit of small C beyond the measured range, we would 
expect the blocked current to recover J^ ocC, as the Don- 
nan potential exceeds the original barrier C/(0). 

Let us switch to the current-voltage characteristics 



shown in the inset of Fig. 3. It is super-linear, because 
the transport is limited by electrostatic barriers shown on 
Fig. 7. On the other hand, this barrier is relatively flat 
and therefore the super-linearity should be weak. This 
agrees with inset of Fig. 3. 

Finally let us concentrate on the sharp dependence of 
DNA capture rate on C (Fig. 2J. It was suggested ^^ 
that DNA loses conformational entropy during translo- 
cation, and therefore the DNA capture rate acquires a 
barrier. However, this mechanism cannot explain the 
observed C dependence, because the ssDNA persistence 
length decreases with C , making the conformation bar- 
rier larger. We notice that the screening cloud of the 
piece of DNA in the channel is squeezed and correspond- 
ing loss of entropy represents a barrier for the DNA cap- 
ture. This barrier decreases with C qualitatively explain- 
ing the growing capture rate with C . 

To summarize, we studied the dependence of the 
blocked ion current during DNA translocation on ion 
concentration. Our experimental results show non- 
monotonic behavior, which cannot be accounted for by 
steric blockade alone. We show that these results may be 
expalined by the inclusion of the electrostatic energy of 
ions inside the pore. 
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